COMPUTING PARAMETRIC RATIONAL GENERATING 
FUNCTIONS WITH A PRIMAL BARVINOK ALGORITHM 

MATTHIAS KOPPE AND SVEN VERDOOLAEGE 

Abstract. Computations with Barvinok's short rational generating 
functions are traditionally being performed in the dual space, to avoid 
the combinatorial complexity of inclusion-exclusion formulas for the in- 
tersecting proper faces of cones. We prove that, on the level of indicator 
functions of polyhedra, there is no need for using inclusion-exclusion for- 
mulas to account for boundary effects: All linear identities in the space 
of indicator functions can be purely expressed using half-open variants 
of the full-dimensional polyhedra in the identity. This gives rise to a 
practically efficient, parametric Barvinok algorithm in the primal space. 



1. Introduction 

We consider a family of polytopes P q = {x£ R d : Ax. < q } parameter- 
ized by a right-hand side vector q E Q C R m , where the set of right-hand 
sides is restricted to some polyhedron Q. For this family of polytopes, we 
define the parametric counting function c: Q — > N by 

c ( q ) = #(p q nz d ). (1) 

Note that this includes vector partition functions c(A) = #{x € N rf : A'x. = 
A} as a special case. It is well-known that the counting function (1) is 
a piecewise quasipolynomial function. We are interested in computing an 
efficient algorithmic representation of the function that allows to efficiently 
evaluate c(q) for any given q. This paper builds on various techniques 
described in the literature, which we review in the following. 

1.1. Barvinok's short rational generating functions. The foundation 
of our method is an algorithmically efficient calculus of rational generating 
functions of the integer points in polyhedra developed by Barvinok [2] ; see 
also [4]. Let P = P q C H d be a rational polyhedron. The generating 
function of P n Z d is defined as the formal Laurent series 
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using the multi-exponent notation z 



a 



Wi = i^ 1 . If P is bounded, gp is 



a Laurent polynomial, which we consider as a rational function gp. If P is 
not bounded but is pointed (i.e., P does not contain a straight line), there is 
a non-empty open subset U C C d such that the series converges absolutely 
and uniformly on every compact subset of U to a rational function gp. If 
P contains a straight line, we set gp = 0. The rational function gp £ 
Q(zi, . . . , Zd) defined in this way is called the rational generating function 

of p n z d . 

By Brion's Theorem [6], the rational generating function of a polyhedron 
is the sum of the rational generating functions of its vertex cones. Thus the 
computation of a rational generating function can be reduced to the case of 
polyhedral cones. Moreover, the mapping P i— ► gp is a valuation: Let [P] 
denote the indicator function of P, i.e., the function 



The valuation property is that any (finite) linear identity X^e/ = 

with £j E Q carries over to a linear identity 'Yli&i £ i9Pi(' L ) = 0- Hence, 
it is possible to use the inclusion-exclusion principle to break a polyhedral 
cone into pieces and to add and subtract the resulting generating functions. 
Indeed, by triangulating the vertex cones, one can reduce the problem to 
the case of simplicial cones. 

By elimination of variables it suffices to consider the case of full-dimen- 
sional simplicial cones, i.e., cones C C R rf generated by d linearly indepen- 
dent ray vectors bi, . . . , £ Z d . The index of such a cone is defined as the 
index of the point lattice generated by bi, . . . , b^ in the standard lattice Z d ; 
we have indC = |det(bi, . . . ,b^)|. Using Barvinok's signed decomposition 
technique, it is possible to write a cone as 



with at most d full-dimensional simplicial cones Cj of lower index in the sum 
over i £ I\ and 0(2 d ) lower-dimensional simplicial cones C% in the sum over 
i £ I-2- The lower-dimensional cones arise due to the inclusion-exclusion 
principle applied to the intersecting faces of the full-dimensional cones. The 
signed decomposition is then recursively applied to the cones Cj, until one 
obtains unimodular (index 1) cones, for which the rational generating func- 
tion can be written down trivially. Since the indexes of the full-dimensional 
cones descend quickly enough at each level of the decomposition, one can 
prove the depth of the decomposition tree is doubly logarithmic in the in- 
dex of the input cone. This gives rise to a polynomiality result in fixed 
dimension: 

Theorem 1 (Barvinok [2]). Let the dimension d be fixed. There exists a 
polynomial-time algorithm for computing the rational generating function of 
a polyhedron P C R d given by rational inequalities. 

Despite the polynomiality result, the algorithm was widely considered to 
be practically inefficient because too many, 0(2 d ), lower-dimensional cones 
had to be created at every level of the decomposition. Later the algorithm 




[C] = ^ [Ci] + £>i PA With £i € ^ 
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was improved by making use of Brion's "polarization trick" , see [6] and [4, 
Remark 4.3]: The computations with rational generating functions are in- 
variant with respect to the contribution of non-pointed cones (cones contain- 
ing a non-trivial linear subspace). The reason is that the rational generating 
function of every non-pointed cone is zero. By operating in the dual space, 
i.e., by computing with the polars of all cones, lower-dimensional cones can 
be safely discarded, because this is equivalent to discarding non-pointed 
cones in the primal space. Thus at each level of the decomposition, only 
at most d cones are created. This dual variant of Barvinok's algorithm has 
efficient implementations in LattE [8, 9, 10] and the library barvinok [19]. 

1.2. Parametric polytopes and generating functions. The vertices of 
a parametric polytope P q = {x € R d : Ax < q}, with q£QC R m are 
affine functions of the parameters q and can be computed as follows. A set 
B of d linearly independent rows of the inequality system Ax < q is called 
a simplex basis. The associated basic solution x(B) is the unique solution 
of the equation Abx = q#. Note that different simplex bases may give rise 
to the same basic solution. A simplex basis (and the corresponding basic 
solution) is called (primal) feasible if Ax(B) < q holds for some q G Q. The 
vertices of P q correspond to the feasible basic solutions and they are said to 
be active on the subset of Q for which the basic solutions are feasible. 

A chamber of the parameterized inequality system Ax < q is an inclusion- 
maximal set of right-hand side vectors q that have the same set of primal 
feasible simplex bases. The chamber complex of P q is the common refine- 
ment of the projections into Q of the n-faces of the polyhedron P = { (x, q) G 
R d x Q : Ax < q}, where n is the dimension of the projection of P onto 
Q [15, 20]. Alternatively, the problem may be translated into a vector parti- 
tion problem, for which the chambers can be computed either directly [1] or 
as the regular triangulations of its Gale transform [12, 17]. However, these 
alternative computations, discussed in more detail in [11, 19], may lead to 
many chambers that do not meet Q and that hence have to be discarded. 

Within each (open) chamber of the chamber complex, the combinatorial 
type of Pq remains the same and Barvinok's algorithm can be applied to 
the vertices active on the chamber [4, Theorem 5.3]. As we will explain 
in more detail in Section 3.1, the result is a parametric rational generating 
function where the parameters only appear in the numerator. In practice, it 
is sufficient to apply Barvinok's algorithm in the closures of the chambers of 
maximal dimension [7, Section 4.2]. On intersections of these closures one 
obtains possibly different representations of the same parametric rational 
generating function. 

Example 2. As a trivial example, consider the one-dimensional parametric 
polytope P q = { x G R 1 : x > 0, 2x < q + 6, x < q}. Its vertices are 0, 
q/2 + 3 and q, active on { g > }, {q > 6} and { q < 6 }, respectively. The 
full-dimensional (open) chambers are { < q < 6 } and { q > 6 } and the 
resulting parametric counting function is 




if < q < 6 
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As in the non-parametric case, P q can be assumed to be full- dimensional 
for all parameter values in the chambers of maximal dimension. Note that a 
reduction to the full-dimensional case may involve a reduction of the param- 
eters to the standard lattice [16, 22]. This parametric version of the dual 
variant of Barvinok's algorithm has also been implemented in barvinok [19] 
and is explained in more detail in [20, 21, 22]. 

1.3. Irrational decompositions and primal algorithms. Recently, Beck 
and Sottile [5] introduced irrational triangulations of polyhedral cones as a 
technique for obtaining simplified proofs for theorems on generating func- 
tions. Let v + C C H d be a full-dimensional affine polyhedral cone; it can 
be triangulated into simplicial full-dimensional cones v + C%. Then there 
exists a vector v € R rf such that 

(v + c) n z d = (v + c) n z d (2) 

and 

0(v + Ci) n Z d = 0, (3) 

that is, the affine cones v + d do not have any integer points in common. 
Thus, without using the inclusion-exclusion principle, one obtains an iden- 
tity on the level of generating functions, 

g v+c (z) = gv +c (z) = ^cfr +c ..(z). (4) 

i 

Koppe [13] considered both irrational triangulations and irrational signed 
decompositions. He constructed a uniform irrational shifting vector v which 
ensures that (3) holds for all cones v+Cj that are created during the course of 
the recursive Barvinok decomposition method. The implementation of this 
method in a version of LattE [14] was the first practically efficient variant 
of Barvinok's algorithm that works in the primal space. 

The benefits of a decomposition in the primal space are twofold. First, 
it allows to effectively use the method of stopped decomposition [13], where 
the recursive decomposition of the cones is stopped before unimodular cones 
are obtained. For certain classes of polyhedra, this technique reduces the 
running time by several orders of magnitude. 

Second, for some classes of polyhedra such as the cross-polytopes, it is 
prohibitively expensive to compute triangulations of the vertex cones in the 
dual space. An all-primal algorithm [13] that computes both triangulations 
and signed decompositions in the primal space is therefore able to handle 
problem instances that cannot be solved with a dual algorithm in reasonable 
time. 

1.4. The contribution of this paper. The irrationalization technique of 
[5, 13] can be viewed as a method of translating an inexact identity (i.e., 
an identity modulo the contribution of lower-dimensional cones) of indicator 
functions of full-dimensional cones, 

Si [vj + C{] = (mod lower-dimensional cones) (5) 
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to an exact identity of rational generating functions, 



(6) 



We remark that this identity is not valid on the level of indicator functions. 
In contrast, in Section 2.1 we provide a general constructive method of 
translating an inexact identity (5) of indicator functions of full-dimensional 
cones to an exact identity of indicator functions of full-dimensional half-open 
cones, 



without increasing the number of summands in the identity. 

This general result gives rise to methods of exact polyhedral subdivision of 
polyhedral cones (Section 2.2) and exact signed decomposition of half-open 
simplicial cones (Section 2.3). 

Since the rational generating function of half-open simplicial cones of 
low index can be written down easily (Section 3.1), we obtain new primal 
variants of Barvinok's algorithm. The new variants have simpler imple- 
mentations than the primal irrational variant [13, Algorithm 5.1] and the 
all-primal irrational variant [13, Algorithm 6.4] because computations with 
large rational numbers can be replaced by simple, combinatorial rules. 

The new variants based on exact decomposition in the primal space are 
particularly useful for parametric problems. The reason is that the method 
of constructing the half-open polyhedral cones only depends on the facet 
normals and is independent from the location of the parametric vertex. In 
contrast, the irrationalization technique needs to shift the parametric vertex 
by a vector s which needs to depend on the parameters. This is of particular 
importance for the case of the irrational all-primal algorithm, where the 
irrational shifting vector s needs to be constructed by solving a parametric 
linear program. 

Moreover, the technique of exact decomposition can also be applied to 
the parameter space Q, obtaining a partition into half-open chambers Qj. 
This gives rise to useful new representations of the parametric generating 
function <?p q (z) (Section 3.2) and the counting function c(q) (Section 3.3). 
We also introduce algorithmic representations of gp q (z) and c(q) that make 
use of half-open activity domains of the parametric vertices. Its benefit is 
that it is of polynomial size and has polynomial evaluation time even when 
the dimension m of the parameter space varies. 

Taking all together, we obtain the first practically efficient parametric 
Barvinok algorithm in the primal space. 

2. Exact triangulations and signed decompositions into 
half-open polyhedra 

2.1. Identities in the algebra of indicator functions, or: Inclusion 
exclusion is not hard for boundary effects. We first show that iden- 
tities of indicator functions of full-dimensional polyhedra modulo lower- 
dimensional polyhedra can be translated to exact identities of indicator 
functions of full- dimensional half-open polyhedra. 




(7) 
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Theorem 3. Let 

5>[Pi] + 5>Fi]=0 (8) 

ieh ieh 

be a (finite) linear identity of indicator functions of closed polyhedra Pi C 
H d , where the polyhedra Pi are full- dimensional for i G I\ and lower- 
dimensional for i G 1%, and where £j G Q. Let each closed polyhedron be 
given as 

Pi = { x : (b* , x) < /3ij for j G J, } . (9) 

Let y € R d &e a vector such that (b* •, y) 7^ /or a// i £ I\\J 1%, j £ Ji- For 
i € Ii, we define the half-open polyhedron 



Pi = x 6 R : (b*,,x) < (3 Lj for j G J, with (b* ,,y) < 0, 

1 , (10) 

(bf^x) < A,j for j G Ji with (b^,y) > j. 

JTien 

J>[P]=0. (11) 

Proof. We will show that (11) holds for an arbitrary x G R rf . To this end, 
fix an arbitrary x G R d . We define 

xa = x + Ay for A G [0, +00). 

Consider the function 



/: [O.+oo) 9Ah ^ ei [Pi])(xA). 



We need to show that /(0) = 0. To this end, we first show that / is constant 
in a neighborhood of 0. 

First, let i G L\ such that x G P%. For j G Ji with (b* - ,y) < 0, we have 
(b^-.x) < Pij, thus (b^,x A ) < /3ij. For j G Ji with (bj^-.y) > 0, we have 

(b*j-,x) < thus (b* - ,x A ) < for A > small enough. Hence, x\ G Pi 
for A > small enough. 

Second, let i G Ji such that x ^ Pj. Then either there exists a j £ 
with (b^y) < and (b^x) > ft,.,-. Then (b^.,x A ) > fa for A > 
small enough. Otherwise, there exists a j £ Ji with (b* -,y) > and 
(b?-,x) > Pij. Then (h*j,x\) > (3ij. Hence, in either case, x\ ^ Pi for 
A > small enough. 

Next we show that / vanishes on some interval (0, Ao). We consider the 
function 

g: [0,+oo) 9Ah (j>i[P] + J>i[P]) (x A ) 

which is constantly zero by (8). Since [Pj](x A ) for i £ I2 vanishes on all but 
finitely many A, we have 



g(X)= (^>[P]Vx A ) 



PRIMAL PARAMETRIC BARVINOK ALGORITHM 



V 



for A from some interval (0, Ai). Also, [Pj](x,\) = [Pj](x^) for some interval 
(0, A2). Hence /(A) = g(A) = for some interval (0, Ao). 

Hence, since / is constant in a neighborhood of 0, it is also zero at A = 0. 
Thus the identity (11) holds for x. □ 

Remark 4. Theorem 3 can be easily generalized to a situation where the 
weights Si are not constants but continuous real-valued functions. In the 
proof, rather than showing that / is constant in a neighborhood of 0, one 
shows that / is continuous at 0. 

2.2. The exact polyhedral subdivision of a closed polyhedral cone. 

For obtaining an exact polyhedral subdivision of a full-dimensional closed 
polyhedral cone C = cone{bi, . . . , b n }, 

[ci = £[£<], 

we apply the above theorem using an arbitrary vector y G int C that avoids 
all facets of the cones Cj, for instance 



y = V(l + 7% 



i=l 



for a suitable 7 > 0. 



2.3. The exact signed decomposition of half-open simplicial cones. 

Let C C R d be a half-open simplicial full-dimensional cone with the double 
description 

C = { x G R d : (b*, x) < for j G J< and (b*, x) < for j G J< } (12) 
C = I X j h j : > for j G J< and Xj > for j G J< } (13) 

where J< U J< = {1, . . . , d}, with the biorthogonality property for the outer 
normal vectors bi- and the ray vectors bj, 

(b*,b l ) = -k J = r 1 i \i =j : (14) 

I U otherwise. 

In the following we introduce a generalization of Barvinok's signed decompo- 
sition [2] to half-open simplicial cones Cj, which will give an exact identity 
of half-open cones. To this end, we first compute the usual signed decom- 
position of the closed cone C = cl C, 

[C] = ^^£j[Ci] (mod lower-dimensional cones) (15) 

i 

using an extra ray w, which has the representation 

d 

w = ^^ajbj where = —(b*, w). (16) 
i=i 

Each of the cones C, is spanned by d vectors from the set {bi, . . . , b^, w}. 
The signs G {±1} are determined according to the location of w, see [2]. 
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An exact identity 

[C] =^2 £ i[Ci] witheG{±l}, 

i 

can now be obtained from (15) as follows. We define cones Cj that are half- 
open counterparts of C{. We only need to determine which of the defining 
inequalities of the cones Cj should be strict. To this end, we first show how 
to construct a vector y that characterizes which defining inequalities of C 
are strict by the means of (10). 

Lemma 5. Let 

* i = l 1 f0rieJ ^> (17) 
[-1 forieJ <} 

and let y G R = int cone{ cribj, . . . , cr^b^} be arbitrary. Then 

J< = {jG{l,...,d}:(b*,y)<0}, 

J< = [j € {l,...,d} : (b*,y) >0}. 

We remark that the construction of such a vector y is not possible for a 
half-open non-simplicial cone in general. 

Proof of Lemma 5. Such a y has the representation 

y = Ajbj - ^2 witn > °- 

ieJ< «eJ< 



Thus 

(b*,y) 



J -Xj for j G J<, 
y+Xj forj€J<, 

which proves the claim. □ 



Now let y G R be an arbitrary vector that is not orthogonal to any of 
the facets of the cones Cj. Then such a vector y can determine which of the 
defining inequalities of the cones Cj are strict. 

In the following, we give a specific construction of such a vector y. To this 
end, let b m be the unique ray of C that is not a ray of C{. Then we denote 
by b^ m the outer normal vector of the unique facet of not incident to w. 
Now consider any facet F of a cone Cj that is incident to w. Since Ci is 
simplicial, there is exactly one ray of Cj, say b/, not incident to F. The 
outer normal vector of the facet is therefore characterized up to scale by the 
indices / and m; thus we denote it by b^ m . See Figure 1 for an example of 
this naming convention. 

Let bo = w. Then, for every outer normal vector b^ m and every ray bj, 
i = 0, . . . , d, we have 



0i;l,m ■= -(b* >m ,bi) < 



> for i = I, 

= fori^Z,m, (18) 
G R for i = m. 
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Now the outer normal vector has the representation 

d 

i=l 

The conditions of (18) determine the outer normal vector b^ m up to scale. 
For the normals bQ m , we can choose 

t>o,m = Oi m h* m . (19) 

For the other facets h* lm , we can choose 

Km = \ a m\ K - signa m • aih* m . (20) 

Now consider 

d 

y = X)°'i(l«il+7 i )bi, (21) 

i=l 

which lies in the cone R for every 7 > 0. We obtain 

(b^ m ,y) = -<J m a m {\a m \ + 7 m ) (22) 

and 

(b*, m ,y) = (bj,y) - signa m • Q;(b^,y) 

= - |a m |<r;(|a;| +7') + signa m • aia m {\a m \ +7") 
= (s\gn(aia m )a m - a{) \oti\ \a m \ 

- 07 |am| 7* + sign(a ; a m )<7 m \a{\ 7 m , (23) 

for / 7^ 0. The right-hand side of (23), as a polynomial in 7, only has finitely 
many roots. Thus there are only finitely many values of 7 for which a scalar 
product (b^* m , y) can vanish for any of the finitely many facet normals b^* m . 
Let 7 > be an arbitrary number for which none of the scalar products 
vanishes. Then the vector y defined by (21) determines which of the defining 
inequalities of the cones Cj should be strict. 

Remark 6. It is possible to construct an a-priori vector y that is suitable 
to determine which definining inequalities are strict for all the cones that 
arise in the hierarchy of triangulations and signed decompositions of a cone 
C = conejbi, . . . ,b n } in Barvinok's algorithm. The construction uses the 
methods from [13]. Let < r € Z and y € and such that the open 
cube y + S tX) (^) is contained in C . (For instance, choose y = Y17=i b« an d 
choose r large enough.) Let D be an upper bound on the determinant of 
any simplicial cone that can arise in a triangulation of C, for instance 

D= (max? =1 llbill 2 ) (24) 



by Lemma 16 of [13]. Let C = max" =1 1 1 1*« 1 1 00 - Usi n g t ne data from Theo- 
rem 11 of [13], 



k 



1 log 2 log 2 D 



M = 2(d-l)\(d k C) d -\ 
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we define 



1/1 1 1 

s = - 



r 



(2M) 1 ' (2M) 2 ''"' (2M) d 



Finally let y = y + s. Then (b*, y) 7^ for any of the facet normals b* that 
can arise in the hierarchy of triangulations and signed decompositions of the 
cone C. 

Remark 7. For performing the exact signed decomposition in a software 
implementation, it is not actually necessary to construct the vector y and 
to evaluate scalar products. In the following, we show that we can devise 
simple, "combinatorial" rules to decide which defining inequalities should 
be strict. To this end, let 7 > in (21) be small enough that none of the 
signs 

o"i,m = - sign(b* m , y) 

given by (23) change if 7 is decreased even more. We can now determine 
ai >m for all possible cases. 

Case 0: a m = 0. The cone would be lower-dimensional in this case, since 
w lies in the space spanned by the ray vectors except b m , and is hence 
discarded. 

Case 1: 1 = 0. From (22), we have 

o"o,m = sign(a m )cr m . 
Case 2: I 7^ 0, a\ = 0, a m 7^ 0. Here we have (b* , y) = — 07 \a m \ j l , thus 

07, m = <?i- 

Case 3: I 7^ 0, aia m > 0. In this case (23) simplifies to 

tf>l m ,y) = (o"m ~ CFl) Wl\ |Om| ~ °l \(Xm\ J 1 + °m \&l\ 1™ ■ (25) 

Case 3 a: ai = a m . Here the first term of (25) cancels, so 



cri, m = -sign(b* ,y) 



1 if I < m, 
— 1 if I > m. 



Case 3b: 07 7^ a m . Here the first term of (25) dominates, so 

c%m = -sign(bf m ,y) = 07. 
Case 4 : I 7^ 0, ceictm < 0. In this case (23) simplifies to 

(bz* m ,y) = -ipm + 07) \a m \ - 01 \a m \-f l - a m \a t \ 7" 1 . (26) 
Case 4 a: a\ = a m . Here the first term of (26) dominates, so 

&l,m = <?l = &m- 

Case 4 b: 07 7^ a m . Here the first term of (26) cancels, so 



a lm = -sign(bf m ,y) 



if / < m, 
if / > m. 
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Figure 1. Signed decomposition of a half-open 3- 
dimensional simplicial cone. Each cone is represented by a 
vertex figure. Closed facets (a = 1) are shown in solid lines, 
while open facets (a = — 1) are shown in broken lines. 



Example 8. Consider (the vertex figure of) the three-dimensional cone on 
the left of Figure 1. The open and closed facets can be described as 

ax = -l <J 2 = l cr 3 = -1, 

while the extra ray w = Yli=i a J°i ls such that 

ati < «2 > 03 > 0. 

For the facets of the cones in the decomposition we have 

00,3 = 03 = — 1 00,2 = (T 2 = 1 <7 0j l = —01 = 1 

4a 1 46 4b 1 

Cl,3 = 0"1 = — 1 01,2 = 01 = — 1 02,1 = 01 = — 1 

3b -i 3b 4a .. 

02,3 = 02 = 1 03,2 = 03 = — 1 0'3,1 = 0'1 = — 1- 

The result is shown on the right of Figure 1. 

Finally we remark that other constructions of y are possible, giving rise to 
different combinatorial rules. For instance, the implementation barvinok [19] 
uses a set of rules that correspond to a modification of (21), where for all 
i G J< the coefficient 7* is replaced by ^ l+d . 

3. Parametric Barvinok algorithm using exact decompositions 

in the primal space 

In the previous section, we have shown how to both triangulate a closed 
polyhedral cone (Section 2.2) and apply Barvinok's decomposition (Sec- 
tion 2.3) in the primal space without introducing (indicator functions of) 
lower-dimensional polytopes. The result is a signed sum of half-open simpli- 
cial cones. The final remaining step in obtaining a generating function for 
a polytope is therefore the computation of the generating function of such 
a cone. 

3.1. The generating function of a half-open simplicial rational cone. 

If v(q)+C is a closed simplicial affine cone where C = { Y^j=i ^j^j : > } 
with bj € Z d , then it is well known [18] that the generating function g v ^) + c 
of v(q) + C is 
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with the fundamental parallelepiped of v(q) + C, 

n = v(q) + {^A i b i :0<A i <l}. 

1 3=1 } 

For a half-open cone v(q) + C given by (13), the same formula holds with 

n = v(q) + { E A i b J : < Aj < 1 for j G J< and < Aj < 1 for j G J< j. 

I j=i J 

To enumerate all points in II n Z rf and compute the numerator of (27), 
we follow the technique of [3, Lemma 5.1], which we adapt for the case of 
half-open cones. 

Lemma 9. Let B be the matrix with the hj as columns and let S be the Smith 
normal form of B, i.e., BV = WS, with V and W unimodular matrices and 
S a diagonal matrix S = diags. Then, 

IT n Z d = { a(k) : kj € Z, < kj < sj }, 

with 

a(k)=v(q)+ {<b;,v(q)-WTc»b i + ^ {{<b*, v(q) - Wk)}} bj 

= wk-J2 L< b l> v (q) - ^ k >J b J - E K b ^ v (q) - wk ) - ^ b ^ 

ieJ< jeJ< 

mi/i {•} i/te (lower) fractional part {x} = x — \x\ and {{•}} the (upper) 
fractional part {{x}} = x — \x — 1] = 1 — {— x}. 

Proof. It is clear that each o:(k) G IlnZ^. To see that all integer points in II 
are exhausted, note that deti? = det5" and that all a(k) are distinct. The 
latter follows from the fact that cc(k) can be written as a(k) = W\t + B~f = 
Wk + WSV -1 ^ for some 7 G Z d . If a(ki) = cc(k2), we must therefore have 
ki = k2 (mod s), i.e., ki = k2. □ 

3.2. Representations of the generating function of a parametric 
polytope. Let Q±, . . . ,Qk Q Q be the chambers of the parameterized in- 
equality system Ax < q of maximal dimension. For all parameters q from 
any given chamber Qi, the parametric polytope P q = {x G R d : Ax < q} 
has the same set of primal feasible simplex bases. Due to amne-linear de- 
pendencies in the set Q of parameters, several primal feasible simplex bases 
can yield the same vertex of the polytope P q on the whole chamber Qi. 
By this mapping we obtain a set Vi of parametric vertices v, (q) for j G V% 
and associated vertex cones Vj(q) + Cj. Let us denote by g v -fq)+c ( z ) the 
parametric generating function of the vertex cone at Vj (q) . 
By Brion's Theorem, we obtain the expression 



E fyw+cM ( 28 ) 



for the generating function of the parametric polytope P q , valid for all pa- 
rameters q G Qi. It turns out [7, Section 4.2] that the formula (28) is also 
valid on the closure cl Qi of the chamber Qi. In this way, we obtain the usual 
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representation of the parametric generating function as a piecewise function 
defined on the whole parameter space Q: 

{EjeV! 0vj(q)+a, ( z ) if q G cl Qi 
: (29) 
Eje v k 9v 3 (q)+a, (z) if q G cl Q k . 

As explained in Section 1.2, this yields possibly different expressions for 
values of q on the intersecting boundaries of two or more chambers. 

We are now interested in a different representation of the parametric 
generating function, 

EjeVi^W+C^z) ifqGQi 

<?P» = I • (30) 

k Ejev fc 9 Vj (q)+c, ( z ) if q G Qfc • 

where the sets Qi form a partition of the parameter space, 

Q = Qi U • • • U Q k with ft n Q- = for i / (31) 

The benefit of representation (30) is that it can be rewritten in the form of 
a closed formula using indicator functions, 

k 

5p» = Z)K?<](q) E 3v,( q )+C,(z)- (32) 

i=i jey s 

Clearly, representations (30) and (32) can be obtained by taking the cham- 
bers of all dimensions, since they form a partition of Q. However, we can 
do better: 

Lemma 10. We can construct representations (30) and (32), where k is 
the number of chambers of Ax < q of maximal dimension. When the di- 
mension d of the polytopes and the dimension m of the parameter space are 
fixed, the construction is possible in polynomial time. 

Proof. Again, we can apply the technique of Theorem 3 to define half-open 
polyhedra Qi that satisfy (31), where y is now an arbitrary vector from 
the relative interior of one of the chambers of maximal dimension. The 
complexity in fixed dimensions m and d follows from the fact that there are 
only polynomially many full-dimensional chambers in this case. □ 

Note that the generating function of a parametric vertex may appear 
multiple times in representation (32) since a vertex Vj(q) may be active 
on more than one chamber. The multiple occurrences can be removed by 
considering the activity regions 

Aj ={q:>lv i (q) < q} 

of individual vertices instead of the chambers. Then, by introducing their 
half-open counterparts Aj constructed by Theorem 3, we obtain another 
representation of the parametric generating function, 

9P^( Z ) = ^[^](q)5 , v J (q)+c* J ( z )- (33) 
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where V = V± U ■ ■ ■ U is the index set of all appearing parametric ver- 
tices. One advantage of this representation is that it can be computed in 
polynomial time, even if the dimension m of the parameter space varies: 

Lemma 11. The representation (33) can be constructed in polynomial time 
when the dimension d of the polytopes is fixed (but the dimension m of the 
parameter space varies). 

Proof. This follows from the above discussion; the number of parametric 
vertices is polynomial when the dimension d of the polytopes is fixed and 
the dimension m of the parameter space varies. □ 

3.3. From the generating function to the counting function. After 
computing the parametric generating function <7p q (z) of P q , an explicit rep- 
resentation of the parametric counting function c(q) = #(-P q H Z d ) can be 
obtained by evaluating the generating function at 1, i.e., c(q) = gp q (l). 
Care needs to be taken in this evaluation since 1 is a pole of each term in 
5-P q ( z )- One typically computes the constant terms of the Laurent expan- 
sions of these rational functions; see [2, 4, 9, 22]. 

Applying this process to (32) and (33), one obtains the counting formulas 

k 

c(q) = J][Qi](q) Yl C v J (q)+C J 

and 

c(q) = J^jKqJCv^qJ+Cj-i 

where Cy^qj+c- i s the sum of the constant terms in the Laurent expansions 
of the terms in ff Vi ( q ) + c 3 ,(z). 

3.4. The resulting algorithms. The complete resulting algorithm, based 
on a chamber decomposition, is shown below. 

Algorithm 12 (Primal parametric Barvinok algorithm). 

Input: full-dimensional parametric polytope P q = {x £ H d : Ax < q}, 

with q£QC R m ; the maximum enumerated cone index £ 

Output: parametric counting function c: Q — > N with c(q) = #(-P q H Z d ) 

(1) Compute the chamber decomposition Q C 2^ of P q and for each 
Qi £ Q of maximal dimension, the corresponding active vertices 
Vi = { Vj(q) }j (see Section 1.2) 

(2) Compute half-open chambers Qi from Qi. 

(3) For each vertex cone Vj(q) + Cj of P q , with Vj(q) € [Jq^q Vi 

(a) Triangulate Cj into half-open full-dimensional simplicial cones 
[Cj] = EjJCj'fc] ( see Section 2.2) 

(b) For each Cjk, apply Barvinok's signed decomposition into half- 
open full-dimensional cones [Cjk] = Yli £ jki[Cjki] °f index at 
most i (see Section 2.3) 

(c) For each Cjki, write down the generating function g v .( a \ + Q. (z) (27) 

of the affine cone v_, (q) + Cjki (see Section 3.1) 

(d) Write down g vM)+C] {v) = J2kT,i £ jkl9 Vj[(i)+ c jkl ^) 
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(4) Write down the generating function gp q (z) (32) of the parametric 
polytope P q (see Section 3.2) 

(5) Specialize the generating function gp q (z) to obtain the counting func- 
tion c(q) = 5p q (l) (see Section 3.3) 

We omit the variation based on activity regions, as it is nearly identical. 
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